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Abstract 

We have recently proposed that the Standard Model Higgs might be responsible for 
generating the cosmological perturbations of the universe by acting as an isocurva- 
ture mode during a de Sitter inflationary stage. In this paper we study the level of 
non-Gaussianity in the cosmological perturbations which are inevitably generated 
due to the non-linearities of the Standard Model Higgs potential. In particular, for 
the current central value of the top mass, wc find that a future detection of non- 
Gaussianity would exclude the detection of tensor modes by the PLANCK satellite. 



1 Introduction 



The CERN collaborations ATLAS [1] and CMS [2] have recently reported the discovery of a 
boson, with mass around (125 — 126) GeV, whose branching ratios are fully consistent with the 
Standard Model (SM) Higgs boson. Motivated by this achievement, some of us have recently 
pointed out [3] that the SM Higgs field might be responsible for generating the cosmological 
perturbations of the universe by acting as an isocurvature mode during a de Sitter inflationary 
stage [4]. 

This idea is based on the fact that during a period of exponential acceleration with Hubble 
rate H all scalar fields with a mass smaller than H are inevitably quantum-mechanically excited 
with a final superhorizon fiat spectrum. The comoving curvature perturbation, which provides 
the initial conditions for the Cosmic Microwave Background (CMB) anisotropics and for the 
Large Scale Structure (LSS) of the universe, may be therefore generated on superhorizon scales 
when the isocurvature perturbation, which is associated to the fluctuations of these light scalar 
fields, is converted into curvature perturbation after (or at the end) of inflation [5-10]. In the 
inhomogeneous decay rate scenario [8] the field responsible for inflation, the inflaton, decays 
perturbatively with a decay rate F. The hot plasma generated by the inflaton decay relics 
thermalizes with a reheating temperature ~ (MpiF)^/^. As a consequence, if F depends on 
some light fields having a flat spectrum during inflation, the decay rate will be characterized 
by large-scale spatial fluctuations, thus leading to a temperature anisotropy, 6Tr/Tr ~ 6T/T. 

In Ref. [3] we have therefore assumed that there was an inflationary period of accelerated 
expansion during the primordial evolutionary stage of the universe and that the isocurvature 
perturbations are generated by the SM Higgs field. This allowed us to play with two independent 
parameters, the SM Higgs mass nih and the Hubble rate H. Our findings show that, the Higgs 
can generate the cosmological perturbation in a wide range of values of the Hubble rate during 
inflation; depending on the values of the SM parameters, e.g. the top-quark mass, and for 
the Higgs masses indicated by the recent experimental discovery, we found H = (10^° — 10^^) 
GeV. On the other side, if the forthcoming data from the Planck satellite will provide hints 
of a 5-mode in the CMB polarization originated from tensor modes, this can be recast into a 
well-defined range of the Higgs mass 

,™.)B— . 128.0GeV ^ 1.3 {'^^^^^^^) GeV + 0.9 (^) GeV, ,1, 

establishing a very interesting correlation between collider and cosmological measurements. 

In this paper we further study the idea that the cosmological perturbations may be ascribed 
to the SM Higgs and focus our attention on another interesting observable characterizing the 
perturbations, the so-called non-Gaussianity (NG) [11]. Primordial NG in the cosmological 
perturbations is a key target for current and future observational probes both in the CMB 
anisotropics and in LSS. The reason is simply stated: detecting some level of NG in the pri- 
mordial fluctuations as well as its shape allows to discriminate among different scenarios for the 
generation of the primordial perturbations [11]. This is because the different mechanisms giving 
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rise to the inflationary perturbations are correlated with specific shapes. Models such as the one 
we investigate in this paper, where the curvature perturbation is sourced by some isocurvature 
perturbation, develop the non-linearities when the perturbation is already on super-Hubble 
scales. The resulting NG is a local function of the Gaussian part and therefore is dubbed 
local. As a consequence, the three-point function in Fourier space gets its major contribution 
from the squeezed configuration where one momentum is negligible compared to the others 
{ki ^ A;2 ~ ks). The squeezed limit of NG leads to strong effects on the clustering of dark 
matter halos as the halo bias becomes strongly scale-dependent [12]. A detection of primordial 
NG in the squeezed shape will allow us to rule out all standard single-field models of inflation 
since they all predict very tiny deviations from Gaussianity [13, 14]. 

In this paper we will answer a very simple question: what is the predicted level of NG in 
the cosmological perturbations if the latter are generated by the SM Higgs? Since we do not 
wish to restrict ourselves to any specific mechanism to convert the isocurvature perturbation 
of the SM Higgs into curvature perturbation, we will compute the unavoidable NG generated 
by the non-linearities of the SM Higgs quartic potential. This choice is justified by the fact 
that the amount of NG depending on the specific mechanism of conversion from isocurvature 
to curvature perturbation may be negligible. We will justify this statement with a couple of 
examples. 

The paper is organized as follows. In section 2 we give a short summary of primordial 
NG and how to calculate it. In section 3 we turn our attention to the NG generated by the 
non-linear interactions of the SM Higgs. In section 4 we present our results and draw our 
conclusions. 

2 Some issues about Non-Gaussianities 

In this paper we will adopt the SN formalism [15] to study the curvature perturbation. This 
formalism makes use of the fact that the local expansion in separate regions on sufficiently large 
scales is identified with the expansion of the unperturbed Friedmann. In general, the curvature 
perturbation is given by 



where 6N must be interpreted as the amount of expansion along the worldline of a comoving 
observer from a spatially flat slice at time t^, to a generic slice at time t{, and p and P are 
the local energy and pressure densities, respectively. If the time tf identifies a hypersurface 
of uniform energy density, the curvature perturbation at some time tf can be expressed as a 
function of the values of the relevant light scalar fields (the inflaton field, the SM Higgs, etc.) 
a\t^:,x) at some time t^: 




(2) 



C(tf , x) = Nia' + -Nija'a-^ + ■■■ 



(/ = 1,--- ,M), 



(3) 
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where Nj and Njj are the first and second derivative, respectively, of the number of e-folds 
with respect to the field . The corresponding two-point correlator of the comoving curvature 
perturbation is 

P^ih) = NjNjPi^, 
{alaQ = {2'n)H{k, + k,)Pl'^ ={27:)H{k, + k,)6''^^. (4) 

Notice that, because of the spatial conformal symmetry of the de Sitter geometry, the two-point 
correlators are diagonal in field space [16]. As we are interested in perturbation generated by 
the SM Higgs, we will take M = 2: one field is infiaton cj) and the other the SM Higgs, A^^ ^ Nh-, 
that is the perturbations are dominated by the SM Higgs. 

The three- and four-point correlators of the comoving curvature perturbation, called also 
the bispectrum and trispectrum respectively, can be also obtained from the expansion (3) 

(Cfc.Cg.C^) = Bl-"-{hMM) + Bl^{k,MM. 

4, h) = NjNjkNl {PifPi^ + 2 permutations) (5) 

and 

T--(^„ h, h, k,) = N.NjN^N.Tlfl^^^ + NuNkN,Nm [Plf Bi^^^^^ + 11 permutations) , 
T^-{h,k2,h,k4) = A^/jAT^lATmAT^ (p£Pi*^P^^ + 11 permutations) 

+NjjkNlNmNn (^Pi^Pl^'P^f + 3 permutations) , (6) 

where kij = iki + kj). To take into account the fact that the fields can be NG at Hubble 
crossing, we have defined their three- and four-point correlators 

= (2-)^^(^i + ^^2 + 4 + ^^4)Ti£4,^. (7) 

We see in full generality that the pieces contributing to the three- and the four-point correlators 
of ( belong to two different groups: the first is proportional to the connected correlators of the 
fields; it originates whenever the fields are intrinsically NG. This is certainly true for the 
SM Higgs fluctuations as the SM potential is not quadratic. We have decided to call it the 
non- universal contribution following Ref. [17]. The second universal group originates from the 
modes of the fluctuations when they are super-Hubble, even if the fields are gaussian. 
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In the presence of NG, one can define the following non-linear parameters^ 

/nl = A ^'^^i^^aCca) ^ (squeezed: /ci < /ca ~ ^3) (8) 

TNL = ^^^1^^, (collapsed: ku^O), (9) 

ki ks ki2 

2rNL + ^5'NL = . / . (^fciCfc2Cfc3Cfc4) (squeezed: ^4 < /ci, ^2, ^3)- (10) 

PS (pSpS +2 permutations 

fc4 \ ki k2 J 

These parameters get both non-universal and universal contributions. For instance, 



/ 



f fnon— un _i_ j-un 

JNL — JnL "t iNL' 

5 5™(fci,fc2,fc3) 



non— un 



NL ^2 [NiN^Y 

^ 5 N'N.jN-^ 
•^^^ 6 [NiN^Y ■ ^ ^ 

In this paper we will focus our attention to the contributions to the NG coming from the non- 
universal pieces: we can compute them exactly once the SM Higgs potential is known. The 
universal pieces are model-dependent, i.e. they depend on the specific mechanism by which 
the isocurvature perturbation is converted into the curvature perturbation. We will assume 
that their contribution to the total NG is subdominant, and motivate this assumption with a 
couple of examples in the following two sub-sections. Of course, in a given well-defined model 
one should explicitly check the validity of this hypothesis. 



2.1 Non-Gaussianities in the modulated decay scenario 

As we mentioned in the Introduction, a specific example where the primordial density pertur- 
bations may be produced just after the end of inflation is the modulated decay scenario when 
the decay rate of the infiaton is a function of the SM Higgs field [8], that is F = F(/;,). Therefore, 
from now one we restrict ourselves to the case in which there is only one relevant fluctuating 
field, the SM Higgs. If we approximate the infiaton reheating by a sudden decay, we may find 
an analytic estimate of the density perturbation. In the case of modulated reheating, the decay 
occurs on a spatial hypersurface with variable local decay rate and hence local Hubble rate 
H = T{h). Before the infiaton decay, the oscillating infiaton field has a pressureless equation 
of state and there is no density perturbation. The perturbed expansion reads 

"'^The prime denotes correlators without the (27r)'^(5'^-'(^,- ki) factors. 
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As the universe is made of radiation after the inflation decay, we can write 

C^«. + l,ug), (13) 

Ehminating SN^ and using the local Friedmann equation p ~ H^, to determine the local density 
in terms of the local decay rate F = T{h), we have at the linear order 

1 6T 1 dlnF 6h ^ 6h 
ol b d\nh h h 

where h is the vacuum expectation value of the Higgs field during inflation. The corresponding 
power spectrum of the curvature perturbation is given by 

By going to higher orders in the fluctuation h, we get (primes indicate differentiation with 
respect to the Higgs field) 

fun 

J] 



NL 





F'T\ 




F'2 ) 




25 A^'" 


50 


54 N'^ 





F"F F'^F'"\ 

Now, suppose that the function T{h) gets its dependence on the SM Higgs from some Yukawa- 
type interaction with Yukawa coupling Y = Yo{l + h/2M), with M some high mass scale. 
If so, T{h) ~ Fo(l + h/2MY and we obtain f^l = 5/2 and g^l = 25/3. Both these non- 
linear parameters lead to a level of NG which is difficult to test observationally and therefore 
subdominant. For further discussion about this particular mechanism, see also Ref. [18]. 



2.2 Non-Gaussianities generated at the end of inflation 

Another possibility is that that the dominant component of the curvature perturbation is 
generated at the transition between inflation and the post-inflationary phase [9]. Let us suppose 
that slow-roll inflation suddenly gives way to radiation domination through a waterfall transition 
in hybrid inflation and that the value of the inflaton 0e at which this happens depends on the SM 
Higgs field, 0e = <Pe{h). If we assume again that the contribution to the curvature perturbation 
from the inflaton fleld is subdominant, we get 

C = Ar^50e = Ar^^£ = ^a. (17) 
dlnh h dliah h 

For concreteness, let us assume that inflation ends on those Hubble regions where X^cfP' + g'^h^ = 
constant. In this case, the corresponding power spectrum of linear perturbations is given by 

p_jLri!^VfiV, (18) 
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where eg = —H/H"^ at the time of the waterfall transition. The universal contribution to the 
three-point correlator is given by (primes again indicate differentiation with respect to the SM 
Higgs) 



5^^^ (iV^</.^ 



/■un _ e^e ' e re ^ _|_ „ — 9,^ 




-^^^'^\l + !L__\^ (19) 



h g'^h 

where we have introduced the slow-roll parameter through the relation A^" = N'^{r]e — 2ee). 
If the power spectrum induced by the SM Higgs field dominates over the one generated by the 
inflaton field P^^ = {l/2ekM^{){H/27ry, that is 

(/ VA^0e)' > (ee/ek) > 1 , (20) 
where is the slow-roll parameter when the wavenumber k is exiting the Hubble radius, then 



/---V2i:V?^. (21) 



This value of NG can be small as long as ^/2e^Mpl < h and eg ^ 1, once Eq. (20) is taken into 
account. Having shown that the universal contribution to NG may be small, we now focus our 
attention on the non-universal contribution, which we can compute exactly once the SM Higgs 
potential is known. 



3 The SM Higgs contribution to the non-Gaussianity 



Let us now consider the SM Higgs field h{x, r) with potential V{h) = X{h)h^/4. If the SM Higgs 
field has a non vanishing vacuum expectation value h, we know that the n-point correlator is 
given by 

{hM)h,ir) . . . h^,^iT)) = -^{0 r dr' [h^,^ir)h^,^ir)...h^,^iT),Vihir'))] o). (22) 



(23) 



Using the mode functions for a massless field in de Sitter 

H 



/2F 



[1-tkr) e'"^, 



one obtains [19] 



2n-4 J 



In{kli ^2; ■ ■ ■ ) ^n); 



(24) 
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where k'f'^ ^h + k2 + ---kn, V^'''^ = d''V{h)/dh!' and 



Tend cJt- r („) 

In{kiMr-- .kn)^2 j —Rey-i{l-ik^T)---{l-iknT)e^^ "J. (25) 



In particular, we have 



I3iki,k2,ks) = -- -7T(7E + iVfcJ . ' 



h{ki,k2,k3,k4) = / ' x2 +^~^ ~ o (7E + A^fcJ - - 3 - (26) 

Here 7e denotes the Euler gamma, while N^^ is the number of e-folds from the time the mode 
kt crosses the Hubble radius to the time Tend when inflation ends. 

How large can the non-universal contributions be? Let us consider the three-point correlator 
in the squeezed limit 

(hh^AY ^-I'^-^N,, P^,P^,^, {k, « k2 ~ A:3). (27) 

The corresponding contribution to the three-point correlator of the comoving curvature pertur- 
bation is 

(C.^.C,^.C.^3)' - Pi^l^ (^1 « k2 - ks) (28) 

leading to a non- universal contribution to the non-linear parameter /nl given by (see Eq. (8)) 



/NL - 3^2^/ - e^Hpy^- 10-3 J 1^ 10-3 J [50 J' 

where wc have defined the quantity V<^ = {k^/27r'^)PS = N''^{H/2T^f ~ 2.3 x IQ-^ and the 

(field-dependent) Higgs mass m\{h) = 3Xh . Analogously we find the four-point correlator to 
be 

^(4)^3 
1 « J.3 J,3 1,3 i,3 

XU rv-^ i\i<2 '"^3 ''^4 

which, in the squeezed limit k^ <^ ki,k2,k3, reduces to 

(KKhM - 6A//^ i6fcVfc3i3fe3 ^3(fei, k2. k,) = = ^^{h^hM 

h 

This corresponds to 

- 27(27r)2 - 1^10-3 j 1^50 ' ■ ^ 

Notice that both /nl and (^nl are predicted to be negative. 



{KKh^y = ^>^H \JJJ^^^^j ,{k,MMMl (30) 



''iKKK)'- (31) 
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Figure 1: The absolute value of f-^i, as a function of H , for nih = 125.5 GeV and three values of 
the top mass: mt = 171.7, 172.4, 173.1 GeV. The curves correspond to values of H below the instabiity 
scale, and they end where H = Amst- The limit |/nl| < 10^ has also been displayed as horizontal band. 



4 Results and conclusions 



To find the amount of NG generated by the SM Higgs we need to compute its potential and 
the solution to the field equation of motion in the expanding background. As in Ref. [3], we 
adopt 2-loop renormalization group (RG) equations for gauge, Higgs-quartic and top-yukawa 
couplings; the pole-mass matching scheme for the Higgs and top masses is taken from the 
Appendix of Ref. [20]. We considered mt = 173.1 ± 0.7 GeV for the top mass [22], and 
as{Mz) = 0.1184 ± 0.0007 for the QCD gauge coupling [23]. For simplicity, we have only 
considered the central value of the QCD coupling a^; in fact, the size of the effect of the la- 
variation of as turns out to be comparable with the higher-order effects that we are neglecting 
(e.g. three-loops effects). The RG-improved effective potential V{h) of the SM Higgs [21] is 
then obtained by solving numerically the RG equations, as a function of input parameters, such 
as the top and Higgs masses. 

The dynamics of the SM Higgs field during infiation is found, after setting the Hubble rate 
during infiation and the initial value of the field, by solving the equation 

'ii + 3Hh + V'{h) = , (33) 

where the dot refers to derivative with respect to t = H^^ In a (conventionally, the scale factor 
is set to 1 at the initial field value). The dynamics has to satisfy the following requirements: 
1) V{h) <^ H'^MpY, which means that the contribution of the Higgs field to the energy density 
during infiation is negligible; 2) \d'^V{h)/dh'^\ ^ H^, which means that the SM Higgs field 
is light enough during infiation. The latter condition must hold for enough e-folds to ensure 
sufficient homogeneity and isotropy of the present observable universe (we consider 60 e-folds 



8 



10'5 

10'" 

> 

O 10" 

10'= 
10" 

124.0 124.5 125.0 125.5 126.0 126.5 127.0 

nih [GeV] 

Figure 2: The green hands indicate the allowed values ofmh,H, where H < Ainst and |/nl| < 100, 
for different top-quark masses (the central value, la and 2a below). 

as a fiducial number). An automatic consequence of tlie condition 2 is tliat tlie special index 
of the perturbations is sufficiently close to 1. 

Next, we scan over the Higgs and top masses and the Hubble rate H, looking for those 
values for which the cosmological perturbations have the correct value, P^^'^ = 4.8 x 10~^, and 
computing the amount of NG generated according to Eqs. (29) and (32). Of course, we have 
to impose that the NG are sufficiently small to meet the current bounds, |/nl| ^ 10^ [24] and 
I^nlI < 10^ [25]. 

The instability scale Ainst of the SM Higgs potential at which the Higgs quartic coupling 
runs negative sets an upper limit on the values of H. A lower limit on H is set by the NG 
constraints. The way to keep /nl below the bound is by having X(h) small and this is achieved 
by taking h as close as possible to Ajnst, where A = 0. The field can roll down starting, at most, 
from the maximum of the potential, and this is our initial condition for the equation of motion 
(33). 

Our results are presented in Figs. 1, 2 and 3 where we also show the current exclusion limit 
on H from CMB data which dictates if < 3 x 10^^ GeV [26] (darker gray band), and the 
testable region by the PLANCK experiment 6.7 x 10^^ GeV <H<3x 10^^ GeV [27] (lighter 
gray band). 

In Figure 1, we plot the values of /nl, as a function of H, for fixed nih = 125.5 GeV and 
different values of the top mass: mt = 171.7, 172.4, 173.1 GeV. The curves correspond to values 
of H below the instability scale, and they end where H = Ainst- The limit |/nl| < 10^ has also 
been displayed as a horizontal band. Notice that, for a given ruh, a sufficently small amount 
of NG is produced as long as the top mass is small enough. This because the instability scale 
increases for smaller rut, and this allows larger values of H to be considered, and thus smaller 
/nl- 
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124.0 124.5 125.0 125.5 126.0 126.5 127.0 



mi, [GeV] 

Figure 3: The orange bands indicate the allowed values of m^^H, where H < Ainst (^nd I^nlI < W^, 
for different top-quark masses (the central value, ±la and ±2a). The green bands correspond to the 
regions allowed by the /nl constraint, as in Fig. 2. 

The bands in Fig. 2 correspond to the upper and lower hmits on H discussed above, for 
different values of the Higgs and top-quark masses. In particular, for the central value of 
TTit = 173.1 GeV, H is limited to be in the range (10^^ — 10^^) GeV, and therefore no tensor 
modes would be seen. Notice that the /nl constraint implies a lower bound on H almost 
independent on mh,mt: H > 10^^ GeV. This can be recovered by imposing /nl = —100 in 
Eq. (29). In the range of masses we considered we typically have h > 10^^ GeV and A > 4x 10"*^. 
This implies a minimum value of if ~ 10^^ GeV. 

In Fig. 3, the corresponding limits from (yfNL are shown (orange bands). We deduce that 
whenever /nl is large, ^^nl is small and therefore the bounds on the bispectrum are stronger 
the the ones coming from the trispectrum. 

One can also obtain a fitting formula relating the Higgs mass, the top mass and the Hubble 
rate on the isocurve /nl = —100, 

. 128.1 GeV ^ 1.3 -^^'^l^"^ ) GeV . 1.3 (^) GeV. (34) 

This formula is vahd in the ranges 124.0 GeV < nih < 127.0 GeV and 6.7 x 10^^ GeV < H < 
3 X 10^^ GeV. If tensor modes will be detected with a given H, Eq. (34) makes a prediction for 
the Higgs and top masses corresponding to a level of NG of /nl = —100. 

In conclusion, we have studied the level of NG in the cosmological perturbations which are 
generated by the non-linearities of the SM Higgs potential under the hypothesis that the SM 
Higgs is responsible for the generation of the curvature perturbation (but not for driving infla- 
tion). In particular, we have assumed that the universal contribution to the NG is negligible, 
supporting this hypothesis with a couple of examples. In this sense, the NG studied in this 
paper is unavoidable (barring potential cancellations with the universal contributions). Under 
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these circumstances we have obtained that: 



• The NG in the four-point correlator (in the squeezed hmit) is always negligible once the 
bounds on the amplitude of the three-point correlator are accounted for. 

• The current constraints on /nl imply a lower bound on the Hubble rate during inflation, 
H > 10^^ GeV, in the range of Higgs mass indicated by the LHC experiments. 

• For the current central value of the top mass {rrit = 173.1 GcV) a future detection of NG 
would exclude the detection of tensor modes through the B-mode of the CMB polariza- 
tion. 
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